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ABSTRACT 


Magsat data from the magnetically quiet days of November 6-6, 1979, were 
used to derive a thirteenth degree and order spherical harmonic geomagnetic 

field model, MGST(3/80). The model utilized both scalar and vector data and 
fit that data with standard deviations of 8, 52, 55 and 97 nT for the scalar 
magnitude, B r , B 0 , and B^, respectively. When compared with earlier models, 
the earth's dipole moment continues to decrease at a rate of about 26 nT/yr. 
Evaluation of earlier models with Magsat data shows that the scalar field 
at the Magsat epoch is best predicted by the POGO(2/72) model but that the 
AWC/75 and IGS/75 are better for predicting vector fields. 


2 


INTRODUCTION 

The NASA spacecraft Magsat (NASA/GSFC, 1979) was successfully launched 
Into a low-altitude, near-polar orbit on October 30, 1979. Intended speci- 
fically to map the magnetic field near the Earth (Langel et al., 1977; 

Langel, 1979), Magsat's instrumentation includes a Cesium vapor scalar mag- 
netometer of accuracy +1 nT (nanotesla) and a fluxgate magnetometer of 
accuracy better than 3 nT rms (Acuna et al., 1978). The scalar magnetometer 
has been successfully utilized to calibrate the vector magnetometer to an 
rms difference in field magnitude of less than 1 nT for the data reported 
herein. Initial attitude determination utilized an Earth horizon scanner and 
a precision sun sensor. The resulting accuracy is estimated to be somewhat 
better than 20 arc-minutes rms, which is equivalent to a component accuracy 
of 290 nT in a 50000 nT field or 175 nT in a 30000 nT field. Later data 
processing (Langel and Berbert, 1979) will utilize on-board star cameras and 
an optical system to measure the attitude of the vector magnetometer relative 
to the star cameras. This is expected to improve the attitude accuracy to 
about 20 arc-seconds rms, with corresponding component accuracies of 3-6 nT. 
The accuracy of orbit determination is better than the 60 m radial and 300 m 
horizontal accuracy required for a 3 nT rms overall accuracy of the scalar 
data (Langel, 1976). For the initial attitude accuracy, orbital error has 
negligible effect on the vector data. 

The anticipated lifetime at launch was five months. However, it now 
appears that orbital decay is considerably less than predicted and we expect 
a seven to seven-and-a-half month lifetime before reentry. 
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Initial orbital parameters for Magsat are: 

Apogee: 561.2 km 

Perigee: 352.4 km 

Inclination: 96.76°. 

Several magnetically quiet days of Magsat data have now become available 
and are suitable for deriving a preliminary field model and for evaluating 
existing geomagnetic field models for the epoch of November 1979. 

THE MAGSAT MODEL 

Of the data currently available to us for analysis, November 5 and 6 have 
the lowest Kp indices. The highest three-hourly Kp value is 1+, which occurred 
only once. Plots of Dst and of the data itself also indicated that these days 
were indeed very quiet. Accordingly, a selection cf 7468 data points from 
these days was used to derive an initial Magsat model whose spherical harmonic 
coefficients are given in Table 1. This included 2414 scalar points whose 
distribution is depicted in Figure 1. Global coverage was excellent; all but 
10 of 648 10°xl0° blocks contained at least one data point. Also included were 

1679 values of (L, 1681 values of B. , and 1694 values of B, , where B . B. and 

r 9 <p r 9 

B are the field components measured in spherical coordinates, for latitudes 
below 50°. Their distribution for those latitudes is similar to that of the 
scalar data shown in Figure 1. Vector data above 50° latitude were not included 
because of non-curl-free fields from field-aligned currents. Such fields are 
transverse to the main field and, as shown by Langel (1974), have little or 
no effect on the amplitude of the scalar measurements. All data are at 
twilight local times because of the orientation of the Magsat orbit. 
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The fitting method was essentially that described by Cain et al. (1967). 

A degree and order thirteen model was chosen because it was justified 
by the density of the data and because studies by Phillips (unpublished, 
personal communication), and Cain (1976) indicate that a thirteenth degree 
and order representation is sufficient to represent the main field of the 
Earth. Any such choice is, of course, somewhat subjective. 

This model fits the selected data with the mean and standard deviations 
shown in the following table. 


COMPONENTS 

MEAN 

DEVIATION (nT) 

STANDARD 
DEVIATION (nT) 

Scalar Magnitude 

0 

8 

B r 

22 

52 

B e 

10 

55 

B * 

84 

97 


The source of the non-zero means for the vector data appears to be a 
bias of 4-7 arc-minutes in the initial attitude solution. The standard devia- 
tion of the scalar data is mainly due to a combination of fields from external 
and crustal sources, neither of which is modeled. The standard deviation of 
the vector data is due to current lack of precision in spacecraft attitude. 

Evaluation of such a model is best done by comparison with independent 
data. Although a global set of independent data is not now available to us, 
comparisons have been made with a set of 44 repeat station measurements 
obtained by the U.S. Geological Survey in 1978 and 1979. The root mean square 
(rms) of residuals of total intensity at these stations was determined to 
be of the order of 170 nT, which may be considered as a good fit to the surface 
data. We have also plotted residuals of several additional days of Magsat 


data relative to this model. These plots Indicate that the model Is a good 
representation of the qulet^time field at the spacecraft, 
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The model was used to compute the field magnitude B, declination D, hori- 
zontal field H and vertical field Z at 400 km altitude, Contour plots are 
shown in Figure 2. 

No secular variation terms are included in the model. However, the 
usefulness of the model could be extended in time by including secular varia- 
tion terms from models such as IGS/75, AWC/75, or IGRF 1975 (whose terms con- 
sist of averages of IGS/75 and AWC/75). 

STRENGTH OF THE CENTERED DIPOLE 

The first approximation of the geomagnetic field is that of the field of 
a dipole located at the Earth's center and inclined to its axis of rotation. 

_3 

The strength of that dipole is measured by MR , where M is the dipole 
moment computed from the first three harmonic coefficients and R is the 
mean radius of the Earth. MR has been decreasing for some years (Harwuod 
and Mai in, 1976; Barraclough et al . , 1978), resulting in speculation that we 
may be observing a reversal of the Earth's field. We have computed MR from 
the new model and compared it to earlier models, as summarized in Table 2 and 
Figure 3. The Magsat results show a continuing decrease of MR"' 5 , consistent 
with earlier results. A comparison of the dipole moment measured by Magsat 
with that measured by P0G0 around 1968 indicates a rate of decrease of 26 nT/ 
year. The slope of the line in Figure 3 is -27 nT/year. These rates are 
somewhat larger than the 10-20 nT/year rates measured in the 1940-1960 time 
frame (Harwood and Malin, 1976). If the higher rates continue, the Earth's 
field would reverse in about 1200 years. 
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EVALUATION OF EARLIER MODELS 

The availability of Magsat data makes It convenient to evaluate the pre- 
dictive capability of those field models currently in common use in the 
scientific community (see also Mead, 1979). For this evaluation 4825 scalar 
points and 6637 vector points at all latitudes were selected from 
November 5 and 6, 1979. The best known and perhaps most widely used of the 
current models is the IGRF (International Geomagnetic Reference Field, Zmuda, 
1971; updated for use from 1975, IAGA, 1976). Other models chosen for evalua- 
tion were IGS/75 (Barraclough et al., 1975), AWC/75 (Peddie and Fabiano, 1976) 
and P0G0(2/72) (Langel et al., 1980). Table 3 shows the mean and standard 
deviation of the selected Magsat data with respect to each of these models. 

The P0G0(£/72) model was derived from satellite scalar field data only. 

Its predictive capability for the scalar field at the Magsat epoch is, as 
might be expected, better than the other models shown in Table 3. For com- 
ponent prediction the AWC/75 and IGS/75 models are very close in capability, 
with perhaps a slight edge to IGS/75. The IGRF 1975 model is significantly 
poorer in its predictions of the vector field at the Magsat epoch. 

SUMMARY 

A global set of scalar and vector data from Magsat was used to construct a 
thirteenth degree and order geomagnetic field model at epoch 1979.85. The 
model fits the Magsat scalar data with a standard deviation of 8 nT , and also 
provides a good fit to 44 recent repeat station measurements. Comparison with 
previous models indicates that the Earth's dipole moment continues *o decrease, 
at a rate of 2f nT/year. Evaluation of earlier models with Magsat data shows 
that the P0G0(?/72) model gave the best predictions of the sea la 1 " field, but 
that the IGS/75 and AWC/75 models were better at predicting the vector 


component': . 
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Table 1 
n 

The MGST(3/80) field model. The mean radius of the Earth 
was assumed to be 6371.2 km. Hean epoch Is 1979.85. 
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Table 3: Residuals of selected MAGSAT data to some published field models. Units are nT. 
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Geomagnetic field magnitude (a) and declination (b) at 400 km 
altitude, as computed from the MGST( 3/80) model. 
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